NUCLEON EDM FROM ATOMIC SYSTEMS AND CONSTRAINTS ON 

SUPERSYMMETRY PARAMETERS 
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The nucleon EDM is shown to be directly related to the EDM of atomic systems. From the 
observed EDM values of the atomic Hg system, the neutron EDM can be extracted, which gives a 
very stringent constraint on the supersymmetry parameters. It is also shown that the measurement 
of Nitrogen and Thallium atomic systems should provide important information on the flavor de- 
pendence of the quark EDM. We perform numerical analyses on the EDM of neutron, proton and 
electron in the minimal supersymmetric standard model with CP- violating phases. We demonstrate 
that the new limit on the neutron EDM extracted from atomic systems excludes a wide parameter 
region of supersymmetry breaking masses above 1 TeV, while the old limit excludes only a small 
mass region below 1 TeV. 

PACS numbers: 12.60. Jv,13.40.Em,11.30.Er,14.20.Dh,21.10.Ky,24.80.+y 



I. INTRODUCTION 

In field theory, it is important to understand the prop- 
erty of the symmetry of the Lagrangian density. In the 
standard model, the time reversal invariance is normally 
kept, but due to the phase of Kobayashi-Maskawa mix- 
ing matrix, there is the T-violation term present in the 
Lagrangian density. The phase shows up as the CP vio- 
lation effects on the — mixing in the kaon decay 
processes Since the CP violation of the Kaon decay 

experiment determines the Kobayashi-Maskawa phase, it 
should be interesting to see what should be further physi- 
cal observables which are predicted from the CP violation 
phenomena or T- violating observables. 

Recent experiments on the CP violation of the — 
mixing in the B-meson decay processes suggest that the 
phase of Kobayashi-Maskawa mixing matrix may not be 
able to describe the observed mixing angle from the -B" — 
B^ mixing. This suggests that there may be some other 
terms present in the Lagrangian density, which violates 
the T-invariance. 

The direct measurement of the T-violation effects can 
be examined if one can measure the electric dipole mo- 
ments (EDM) of the isolated system. The finite number 
of the neutron EDM should exhibit most clearly the di- 
rect evidence of the T-violation effects in the Lagrangian 
density. Until now, the upper limit of the neutron EDM 
dn from the direct measurements is around 

d„ ~ (1.9±5.4) X 10"^^ e-cm (1.1a) 

dn ^ (2.6 ± 4.0 ± 1.6) X 10"^^ e • cm. (1.16) 

This EDM is still not finite, but gives a strong constraint 
on the supersymmetry model parameters. 



There are also EDM measurements of atomic systems 
such as i29xe H and ^^^Hg H 0. However, it is 
not very easy to extract the individual particles EDM 
from the atomic EDM measurements since there is the 
Schiff theorem The Schiff theorem states that the 

EDM of the individual particles cannot be measured if 
the constituents are interacting with nonrelativistic static 
Coulomb force. The EDM of the individual particles 
should show u p w hen they are interacting with relativis- 
tic kinematics |llllT^IT3l| or with strong interactions [T^ . 
The EDM due to the former case is found in the heavy 
atomic system while the EDM of the latter case arises 
from the finite nuclear size effects. 

Recent careful studies clarify that the EDM arising 
from the nuclear finite size effects is directly related to 
the neutron EDM dn as 



dxc — 1-Gd„ 



-2.Mn 



(1.2a) 
(1.26) 



where c?xo and dng denote the EDMs of Xe and Hg atomic 
systems. In particular, the EDM of ^^^Hg atomic system 
is observed with high accuracy and from the observed 
vaule of c?Hg, we can deduce the neutron EDM as 



dn ^ (0.37 ± 0.17 ± 0.14) X 10^^* e • cm 



(1.3) 
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This number is three orders of magnitude smaller than 
the direct measurement of neutron EDM in eq.(l.l). 
Therefore, it should be worthwhile to examine how this 
number of the neutron EDM can give constraints on the 
supersymmetry model parameters. 

Supersymmetry (SUSY) is one of the most promising 
candidates for new physics beyond the standard model 
|l6l |. In particular, supersymmetric models contain new 
CP-violating phases |17j in addition to the phase of the 
Cabbibo-Kobayashi-Maskawa mixing matrix, and the ef- 
fect of these phases may explain the discrepancy be- 
tween the experimental result and the standard model 
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prediction for CP-asymmetry in the B-meson decays. In 
general, the new CP-violating phases in supersymmetric 
models induce a non- vanishing EDM for quarks and lep- 
tons. In fact, they give a large contribution to the EDMs 
at the one-loop level. On the other hand, it is known that 
the standard model contributions to the EDMs cancel out 
up to the two-loop level, and the leading three-loop con- 
tribution is much smaller than the current experimental 
upper bound. This implies that a large EDM is not only 
a signal of CP-violation but also an indirect evidence of 
SUSY. 

In the present work, we calculate the EDM of neu- 
tron, proton and electron in the minimal supersymmetric 
standard model (MSSM) with CP-violating phases. We 
shall show that the neutron EDM limit in eq.(1.3) gives 
a very severe constraint on the parameters in the super- 
symmetric model. For a sizable value of the CP-violating 
phases, the magnitude of soft breaking mass parameters 
typically has to be larger than several TeV to satisfy the 
EDM constraint. 

Further, we show, for the first time, that the proton 
EDM can be extracted from the Nitrogen and Thallium 
atomic EDM. This is simply because, for the ^^N case, 
it has the ground state with spin 1/2, and the state is 
described by the proton hole state as |7r(lpi : ^ ) to a 
good accuracy. Also, for the 205x£ case, since it is a single 
proton hole state, it can be expressed as |7r(3si : ^'^) 
to a good approximation even though the neutrons are 
not in the magic but in the two hole state. 

From the experimental constraint, the nucleus should 
have the spin 1/2 such that the quadrupole field should 
not have any effects on the EDM measurements. In this 
respect, ^^N and '^^^T£ atomic systems must be good 
cases for measuring the nuclear EDM. The EDM of the 
Nitrogen atomic system disN and the ^os^p^ atomic sys- 
tem are related to the proton EDM dp, and the cal- 
culated numbers become 

(ii5N =i -O.lGdp (1.4a) 

dTi ~ 4.8dp. (1.46) 

Recently, Regan et.al have measured the atomic EDM 
of by the atomic beam magnetic resonance method 

[T^ . From the measurement of the atomic EDM of ^^^T^, 
we can extract the proton EDM dp 

dp ~ -(0.83 ±0.90) X 10-25 e- cm (1.5) 

which is quite a severe constraint on the proton EDM dp. 
Here, we have not included the electron EDM contribu- 
tion to the atomic EDM of "^^^Ti, and the electron EDM 
contribution will be discussed in a qualitative manner in 
section 4. 

On the other hand, the EDM for the Nitrogen case is 
smaller than the ones in heavier nuclei, but we believe 
that the Nitrogen EDM should be also measured in fu- 
ture. 



This paper is organized as follows. In the next sec- 
tion, we briefly discuss the new mechanism to obtain the 
EDM of the individual particles in nucleus, like neutron 
EDM or proton EDM from atomic systems. Then, in 
section 3, we evaluate the EDM of neutron, proton and 
electron in the minimal supersymmetric standard model 
with CP-violating phases, and find various constraints on 
the supersymmetric parameters from the observed neu- 
tron EDM value. In section 4, we present a qualitative 
discussion on the contributions of the electron and the 
proton EDMs to the atomic EDM value. Section 5 sum- 
marizes what we clarify the nuclear EDM and the related 
parameters which appear in the supersymmetry models. 



II. NUCLEON EDM FROM ATOMIC SYSTEMS 

The electric dipole moment (EDM) of the isolated sys- 
tem can present a direct evidence of the T-violation in- 
teraction if it is finite. Since the measurement is carried 
out for the total atomic system, the EDM of the indi- 
vidual particles (electron EDM or nucleon EDM) should 
be extracted from the atomic EDM measurement. How- 
ever, this extraction is a nontrivial task since there is 
the Schiff theorem. The Schiff theorem states that the 
EDM of the individual particles cannot be measured if 
the neutral system is interacting with nonrelativistic elec- 
trostatic forces. 

The EDM of electron becomes important when the sys- 
tem has a relativistic effect since the Schiff theorem is 
not applicable to the relativistic case. Indeed for heavier 
atoms, the relativistic effect becomes important since the 
relativistic correction is seen as (Za)'^. For the electron 
EDM, there are many calculations in heavy atoms. 

If the interactions are due to nuclear forces, then there 
is no effect from the Schiff theorem. Since the nucleus 
is always found in atomic system, the nuclear EDM can 
be seen as the finite size effects in the total atomic sys- 
tem. Recent calculations on the nuclear EDM show that 
there is an appreciably large effect of the nuclear EDM to 
the total atomic EDM after the effects due to the Schiff 
theorem are taken into account. The results obtained in 
for Xe and Hg atomic EDM are summarized in eqs. 
(1.2). 

Now, we briefly describe the finite size effects of the 
atomic EDM since the detailed calculation is found in 
15] . In particular, the effects due to the Schiff theorem 
are well treated and explained in j^j, and therefore we 
do not repeat it for this part in this paper. 



A. Hamiltonian of atomic systems 

We first write the Hamiltonian of the total atomic and 
nuclear systems for Nitrogen. 

The unperturbed Hamiltonian Hq of the Nitrogen sys- 
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tem can be written 



Pi 

2m 



IE 



2 — ' iTj — r, i 



(2.1) 

where r^, pj denote the coordinate and the momentum 
of the (;lc(;tron while Ri, P, denote the nuclear variable 
and momentum, respectively. 

On the other hand, the perturbed Hamiltonian coming 
from the EDM is written as 

V" V" • (ri - R-i) , V" V" • (r^ - r^) 



i=l jT^i 



- r |3 



1=1 J = l ' ' 



|Ri-R,f 



- ^ • Ee,t - X] • Eext + e ^(ri - Ri) • Ee^t (2.2) 

i=l i=l i=l 

where the summation over Z in nucleus means that it 
should be taken over protons. The EDM of the nucleon 
can be expressed in terms of the nucleon isospin as 

= ^ [(1 + Odp'^' + (1 - <)dncr^] . (2.3) 



B. Nuclear EDM from nuclear excitation 

Here, wc evaluate the finite size effects on the second 
order EDM energy in the nucleus. The perturbed EDM 
Hamiltonian for nuclear parts -fff^^ can be written as 



(l + Tn(R,-R,-) 



|R. 



R.f 



(2.5) 



where and En denote the ground state energy and 
the excitation energy of the nuclear states, respectively. 
The electron states are kept in the ground state through- 
out the calculation, and therefore it is not written here. 
Also, it should be noted that we made use of the relation 

A 

Ri = since we set the center mass coordinate to the 

i=l 

nuclear center. 



C. Nitrogen atomic system 

Now, we calculate the ^^N case in which we assume a 
simple single particle shell model state \tt{1p^)~^). This 
should be rather a good description for ^^N. Further, the 
atomic states stay in the ground state, and therefore the 
electron wave functions are not written here. Thus, we 
write the nuclear wave function as 



I* 



N 



(2.6) 



In this case, the intermediate states \n) that contribute 
to the second order EDM energy [eq.(2.6)] for 7r(lpi ) are 
restricted to 



\n) = |7r(lsi)- 



k(lpl)-^7r(2sl):i^). 

(2.7) 

With eqs. (2.5-7), we carry out numerical calculations of 
the second order EDM energy ^E\^ and obtain the re- 
lation between the Nitrogen EDM dis^ and the proton 
EDM dp as 



5N 



-o.iedp. 



(1.4a) 



Even though the factor in front of dp is not very large 
number, we believe that the EDM of Nitrogen can be 
well measured experimentally in future. 



H. 



in) 
edm 



(2.4) 

Now, we consider the second order EDM energy due to 
the intermediate nuclear excitations, keeping the atomic 
state in the ground state. This process arises from the 
finite nuclear size effects in the EDM Hamiltonian. The 
second order EDM energy can be written as 



AEfJ 



E 



En — En 



i=l 



D. Thallium atomic system 

Now, we calculate the ^^^T£ case in which we assume a 

simple single particle shell model state |7r(3sj.)^^). This 
should be reasonable for the 205rp£ case. Therefore, the 
nuclear wave function for 205^^^ ^^-^ written as 



\^^) = \n{3s.)-':- ). 



(2.8) 



In this case, the intermediate states \n) that contribute 
to the second order EDM energy [eq.(2.6)] for |7r(3si )~^) 
are restricted to 



^ 1 

:{n\Y,j[{l + rndp<T' + {l-rndnCT[ 



|7r(2pi)- 



^ ' '7r(3si) ^7r(3pl) : ^ ). 



), I 



(2.9) 
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With eqs. (2.8-9), we carry out numerical calculations of 

(2) 

the second order EDM energy ^EjJ and obtain the re- 
lation between the Thallium EDM dTi and the proton 
EDM dp as 

dre ^ 4.8dp. (1.46) 

Clearly, the EDM of Thallium is quite large, and there- 
fore there must be a good chance of observing it. Indeed, 
there is a measurement of T£F compound system, and 
Cho el al. ^| obtained the EDM of d^iF as 

dreF ^ (-1-7 ± 2.9) x 10"^^ e • cm. (2.10) 

Making use of the relation between Thallium EDM and 
proton EDM [eq.(1.4b)], we obtain the proton EDM as 

dp ~ (-0.35 ± 0.60) X 10"^^ e • cm (2.11) 

where we ignored the contribution of EDM from the Flu- 
oride atomic system. 

Further, Regan et al. have recently measured the 
atomic EDM of '^'^^T£ by the atomic beam magnetic res- 
onance method pl |. They extracted the electron EDM 
from their measurement, and obtained the constraint of 
electron EDM de- Here, we estimate the proton EDM 
from their measurements, and the atomic EDM of ^os^^ 
is found to be 

dre ^ -(0.4 ± 0.43) x lO'^'' e ■ cm. (2.12) 

Therefore, we can extract the proton EDM using 
eq.(1.4b) 

dp ~ -(0.83 ±0.90) X 10^2^ e- cm (2.13) 

which gives stronger constraints than eq.(2.11) of Cho 
et.al almost by two orders of magnitude. This proton 
value of dp is just comparable to the neutron EDM d„ 
from the direct measurement. 

As mentioned above, the effects due to the electron 
EDM de in heavy atoms should be large. In fact the con- 
tribution of the electron EDM to the Thallium case is 
evaluated in |23|, and they found a very large enhance- 
ment factor. In this paper, however, we have not included 
the electron EDM. Instead, we will discuss the electron 
EDM contribution to the atomic EDM in a qualitative 
fashion in section 4. 



III. CONSTRAINTS FROM NEUTRON EDM 
ON SUPERSYMMETRY MODEL 

In this section, we discuss the constraints from the neu- 
tron EDM in the MSSM. In the present analysis, we work 
in a framework of a phenomenological MSSM. The rel- 
evant parameters in this model are defined as follows. 
SUSY breaking mass terms for the squarks and the slep- 
tons contain a common scalar mass parameter mg as 

-ml{\q\' + \u\^ + \d\^ + \l\' + \i\'), (3.1) 



where q and £ denote the SU(2) doublet fields for the 
left-handed squarks and sleptons, respectively, while u, d 
and e denote the SU(2) singlet fields for the right-handed 
scalar up-quark, the scalar down-quark and the scalar 
electron, respectively. Gaugino mass terms for the bino 
B, the wino (a = 1,2,3) and the gluino {a = 
1, • • • , 8) are given by 

- ^{MiBB + NhW'^W'' + A'hg^'g"') + h.c, (3.2) 

where Mi, M2 and M3 denote a mass parameter for 
U{1)y, SU{2)l and SU{3)c gaugino, respectively. We 
assume the GUT relation for the gaugino masses Mi = 
f tan2 9wM2 and M3 = M2. 

The SUSY conserving part of the Lagrangian can be 
written by the superpotential 

W = Y^U'H2Q + YdD^HiQ + Y,E^HiL + 

(3.3) 

where Y^, Yd and Ye are Yukawa couplings for the up- 
quark, the down-quark and the electron, respectively. 
The coefficient /i is a Higgs mixing mass parameter. The 
symbols with a hat denote chiral superfields with self- 
evident notations. In ea. H3.3|l . we neglect generation mix- 
ings, since we discuss generation conserving effects. The 
scalar trilinear couplings are parameterized by the uni- 
versal mass parameter A as 

A{YuU*H2q + Ydd*Hiq + Yee*Hi£). (3.4) 

Furthermore, the Higgs sector in this model includes 
the ratio of the vacuum expectation values of the two 
neutral Higgs fields, tan/3 = {H2) / {Hi) , and a mass pa- 
rameter niA for the pseudoscalar neutral Higgs field. 

In summary, the Lagrangian is parameterized by the 
following parameters 

tan/3, niA, M2, ttiq, /i, A. (3-5) 

Among these, M2^ /i and A can have a non- vanishing 
complex phase. We take a phase convention that M2 
is real and positive. In order to examine the effects of 
CP-violation, we define two complex phases (0^, 9a) by 

At = |/^|exp(i0^), A = |^|exp(z0^). (3.6) 

The EDMs of quarks and leptons are defined as a co- 
efficient in the electric dipole operator 

C = -'-di!a^,^,JF^r (3.7) 

The EDM d/ of the fermion / {— u, d, e) can be ob- 
tained by evaluating relevant Feynman diagrams. In the 
case of the MSSM, three kinds of diagrams contribute at 
the one-loop level: (i) chargino exchange contribution, 
(ii) neutralino exchange contribution and (iii) gluino ex- 
change contribution. Note that the gluino contribution 



5 



does not exist for the electron EDM de- The expression 
oi df for each fermion can be found in Ref. [2ll| . 

To obtain the EDM of the neutron and the proton, we 
employ the non-relativistic quark model as 

dn = ^(4dd-d«), dp = ^(Mu-dd). (3.8) 

It should be noted that a quark chromomagnetic dipole 
operator and a pure giuonic dipole operator may give 
comparable contributions to the electric dipole operator 
(|3.7|) J2II. Inclusion of these operators may cause can- 
cellations to suppress the neutron EDM. However, this 
kind of cancellations occurs in relatively small portions 
of the whole parameter space. Therefore, we neglect the 
effects of these operators for simplicity. Also, we focus 
on a nuclear finite size effects, and neglect contributions 
from the chromomagnetic dipole operator for the strange 
quark in this analysis. 

In the following, we present the results of our nu- 
merical analysis. In all the figures, we always choose 
tan/3 = 10 and niA = 700 GeV as reference values. For 
large tan/3, the EDMs are typically enhanced so that the 
constraint on the mass parameters becomes more strin- 
gent. The EDMs are not very sensitive to m^, since the 
Higgs masses are not directly related with the relevant 
diagrams. 

Let us begin with the effect of 6a assuming that fi 
is real and positive. Variations of the EDMs in the 
(M2,mo) plane are displayed in Fig. ^ for tan/3 = 10, 
niA = 700 GeV, \fi\ = 500 GeV, |^| = 500 GeV, 9^ = 0, 
9a = 7r/6. The result for the neutron EDM is shown 
in the upper left window. If we assume previous experi- 
mental limit dn ^ 10~^^ using the neutron itself, only the 
small mass region mo, M2 ^ ITeV would be excluded. 
However, with the new limit 

0.06 X lO^^^e • cm < d„ < 0.68 x lO^^^c • cm (3.9) 

derived from eq.(1.3), it is found that the large region 
below the solid line is excluded by the neutron EDM. 

In the upper right window of the same figure, we plot 
the contours of the proton EDM for the same choice of 
the parameters. In this case, dp is negative in the whole 
region. As for the magnitude, \dp\ is of the same order 
of magnitude as dn- This can be seen more clearly in the 
contour plot of the ratio dp/dn (lower left window). In 
the region M2 ^ mo, \dp\ is smaller than dn, while it is 
nearly equal to dn for M2 ^ toq- 

For comparison, we also plot the contours of d^ in 
the lower right window. The current experimental upper 
bound for the electron EDM is \de\ < 1.6 x 10^^''e ■ cm 
[Till . We can see that the constraint from dn is much 
severer than that from d^. 

The effect of varing 9a can be seen in the (mo,^yi) plane 
in Fig. H for tan/3 = 10, ttia = 700 GeV, M2 = 1 TeV, 
= 500 GeV, \A\ = 500 GeV, 9^ = 0. For < 9a < tt, 
dn is positive, and only small regions outside the solid 
lines are allowed by the neutron EDM constraint (|3.9|l . In 



particular, for 0.27r ^ 9a ^ O.Sn, the common sfermion 
mass mo must be larger than 10 TeV to satisfy (|3.9|) for 
this choice of the parameters. For tt < 9a < 2tt, dn is 
negative. 

Finally, we discuss the effect of 0^. The phase de- 
pendence of the neutron EDM for fixed mo and M2 is 
shown in Fig. |31 The relevant supersymmetric parame- 
ters for the upper window are chosen as tan /3 = 10, = 
700 GeV, M2 = 2 TeV, mo = 5 TeV, = 500 GeV, 
1^1 = 500 GeV and 9a = 0, while those for the lower win- 
dow are the same except 6*^ = 0. The range < < tt 
in the upper window and the range tt < 0^ < 27r in the 
lower window corresponds to c?„ < 0. In the region be- 
tween the two dashed lines, the EDM constraint H3.9|) 
is satisfied. From the two plots, we see that the CP- 
violating phases 0^ and 9a are severely bound to zero 
even in a multi-TeV region of SUSY breaking mass pa- 
rameters. Relatively speaking, 9a is less constrained than 
9f^. Therefore, if we take sizable 6*^ rather than 9a, the 
constraints on mo and M2 become severer than those in 
figures ^ a-nd |21 



IV. DISCUSSIONS ON de AND dp 

The expected magnitudes of the electron EDM de and 
the proton EDM dp depend on the parameters of the 
SUSY models. It should be rather difficult to reliably 
predict the relative magnitude between de and dp. 

Therefore, the contributions from de and dp to the 
atomic EDM values should be estimated carefully. In 
this paper, we do not include any contributions of the 
electron EDM, but instead we will discuss the electron 
EDM contribution to the atomic EDM in a qualitative 
fashion. 

Below, we evaluate the contributions from de and dp 
to the Nitrogen and Thallium atomic EDMs. 

A. Electron EDM in atomic system 

The electron EDM contributions to the atomic EDM 
come from the relativistic effects of the electron wave 
functions. This estimation is carefully done in and 
the atomic EDM ^1*°™ can be written as 

^atom _ CiiZaf (1 + 2(Za)2) de. (4.1) 

If we take Ci ^ —100, then we obtain for Nitrogen and 
T£ cases, 

dN ^ -0.34 (4.2a) 

dre ^ -604 (4.25) 

where the sign in front of the number is not determined 
from eq.(4.1) since it depends on each atomic state. 
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The elaborate calculation for T£ case shows ~ -(0.2 - 1.8) X IQ-^^ e • cm. (4.66) 

dre ~ -5854. (4.2c) 



The order of magnitude enhancement may come from the 
atomic state configuration. 



On the other hand, we obtain the atomic EDM for Ni- 
trogen and T£ cases from the proton EDM contribution 
as 



B. Proton EDM in atomic system 

As we discussed in section 2, the proton EDM mainly 
comes from the nuclear excitations. The proton EDM 
contribution to the atomic EDM can be written as 

dz OmZdp. (4.3) 

In this case, we obtain for Nitrogen and T£ cases, 

dN =i ~0.2dp (4.4a) 



~ -0.4 X 10-29 e • cm. (4.6c) 
d% ~ -1.0 X 10-^^ e • cm. (4.6d) 



Therefore, the EDM contributions from the electron (rel- 
ativistic effects) and the proton (nuclear finite size ef- 
fects) for the T£ case may well be comparable. But the 
atomic EDM for the Nitrogen case is mainly determined 
from the proton EDM. 



dre ~ 2.4dp. (4.46) 

On the other hand, as we saw in section 2, the elaborate 
calculations show 



I AT 



-O.Wdp (4.4c) 



V. CONCLUSIONS 



dre ~ 4.8dp. (4.4d) 

As can be seen from eqs.(4.4), the qualitative formula 
of eq.(4.3) gives sufficiently rehable estimations of the 
nuclear EDM. 



C. de and dp from SUSY models 

Now, as we saw in the previous section, the relation 
between the de and dp EDMs from the SUSY model cal- 
culations are shown in fig. 1. If we choose the values of 
M2 and mo to be 

M2 ~ 9 TeV, mo ~ 9 TeV 

then we obtain 

de ~ -2.8 X 10"^^ e • cm (4.5a) 



We have presented the proton and neutron EDMs 
which are extracted from the atomic EDM measure- 
ments. In particular, we show that the proton EDM 
can be obtained if one can measure the EDM of Nitro- 
gen atomic systems. The EDM measurement of Nitrogen 
atomic system enables to discuss the flavor dependence 
of the quark EDM which should give some important in- 
formation on the origin of the T-violation in field theory 
models. 

Also, we have studied the EDM of neutron, proton and 
electron in the minimal supersymmetric standard model. 
In the presence of supersymmetric CP- violation, the new 
limit on the neutron EDM extracted from atomic sys- 
tems excludes a wide parameter region of SUSY breaking 
masses above 1 TeV, while the old limit excludes only a 
small mass region below 1 TeV. Thus the observed neu- 
tron EDM is found to give a very stringent test on the 
supersymmetric model. 



dp ~ -1.9 X 10"^^ e • cm. (4.56) 

In this case, the atomic EDM for Nitrogen and T£ cases 
from the electron EDM become 



-1.0 X 10"^^ e • cm (4.6a) 
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FIG. 1: Contour plots of d„ (upper left window), dp (upper 
right window), dp/dn (lower left window) and de (lower right 
window) in the (M2,mo) plane for tan /3 = 10, niA = 700 GeV, 
ImI = 500 GeV, \A\ = 500 GeV, 6^ =0,9a= 7r/6. 
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FIG. 2: Contours of dn in the (mo, 0a) plane for tan/3 = 10, 
m.A = 700 GeV, Ma = ITeV, |//| = 500 GeV, |^| = 500 GeV, 

0M =0. 
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FIG. 3: The neutron EDM d„ vs Of^ (upper window) and 
Oa (lower window). Relevant supersymmetric parameters for 
the upper window are chosen as tan /3 = 10, mA = 700 GeV, 
Ma = 2 TeV, mo = 5TcV, \n\ = 500 GoV, \A\ = 500 GcV and 
9a = 0, while those for the lower window are the same except 
= 0. The range < < tt in the upper window and 
TT < < 27r in the lower window corresponds to dn < 0. In 
the region between the two dashed lines, the EDM constraint 
0.06 X 10"^*e • cm < d„ < 0.68 x 10"^®e • cm is satisfied. 



